Ab initio Evidence for Strong Correlation Associated with Mott Proximity 

in Iron-based Superconductors 
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We predict that iron-based superconductors discovered near d*" configuration (5 Fe 3ci orbitals filled 
by 6 electrons) is located on the foot of an unexpectedly large dome of correlated electron matter 
centered at the Mott insulator at (namely, half filling). This is based on the many- variable 
variational Monte-Carlo results for ab initio low-energy models derived by the downfolding. The d^ 
Mott proximity extends to subsequent emergence of incoherent metals, orbital differentiations due 
to the Mott physics and Hund's-rule coupling, followed by antiferromagnetic quantum criticality, in 
quantitative accordance with available experiments. 

PACS numbers: 74.70.Xa, 71.27.-|-a, 74.20.Pq ,71.30.+h 



Introduction. — Iron-based superconductors discovered 
in 2008 P soon proved to encompass different families of 
pnictides such as LaFeAsO and BaFe2As2 and chalco- 
genides such as FeSe and FeTe, where the superconduc- 
tivity was discovered mostly under electron or hole dop- 
ing as in the cuprate superconductors The back- 
bone lattices are commonly built from a stacking of iron 
square-lattice layers. The band structures consisting 
mainly of the Fe 2>d orbitals located near the Fermi level 
are also similar among the families Q. However, physi- 
cal properties strongly depend on the families and have 
a diversity in magnetic, transport and superconducting 
properties 

The antiferromagnetic (AF) order is, in most cases, 
found close to the superconductivity, while the ordered 
magnetic moment has a diversity from zero up to > 
2/^B@43 whose origin is not entirely clear. Another 
diversity is found in the coherence of the metallic car- 
rier. The incoherent ( "bad metallic") conduction with 
small Drude weight l9Mlll. enhanced mass in deHaas-van 
Alphen measurement [l2|jStrongly renormalized quasi- 
particle peak in ARPES[13|, unconventional AF fluctu- 



ations [IJI and a typical Mott-Hubbard splitting (emer- 
gence of the lower Hubbard band) [lH have been claimed 
in several compounds, while the correlation effects arc 
less clear in other cases [l6|. Such a rich diversity is char- 
acteristic in contrast to the cuprate superconductors, in 
which strong correlation effects are common and univer- 
sal. In this letter, we first show that the diversity emerges 
indeed from the variation in the electron correlation. 

The strength of electron correlation in the present ma- 
terials is presently under hot debate partly because of 
this diversity. Firm and accurate calculations with full 
account of not only dynamical but also spatial fluctuation 
effects are desired. However, ab initio analyses with full 
account of spatial correlation effects require demanding 
calculations and to our knowledge are very few. 

In this letter, we employ an ab initio method [Till by 
a hybrid combination of density functional theory with 



accurate solvers for the downfolded effective low-energy 
model. Here, we employ the multi-variable variational 
Monte Carlo (mVMC) methodfH, ^ for the solver. 
See Supplemental Materials S.l and S.3 for details of 
the whole methods. The method enables us to examine 
strong correlation effects in an ah initio way fully with 
dynamical as well as spatial fluctuation effects. We find 
that the real iron-based superconductors are on the foot 
of a large-dome structure centered around the Mott 
insulator, whose proximity effects generate various cor- 
relation phenomena of the real superconductors. This 
sheds new light on the understanding of the electron cor- 
relation effects in the iron-based materials. 

In general, to vary and control electronic properties 
in correlated electron systems, two important routes 
are known[l7[. One is the bandwidth (or equivalently 
effective-Coulomb-interaction) control and the other is 
the filling control. The former directly controls the ratio 
between the kinetic and interaction energies of electrons, 
while the latter tunes the distance from the "commen- 
surate filling" (a simple fractional number of band filling 
enhances the electron correlation as in the Mott insula- 
tor) . We show that both play key roles in the iron-based 
superconductors . 

Model. — The derived ah initio parameters for the 
present effective models by the downfolding procedure^, 
l28l | of the Fe Sd orbitals consist of the transfer integrals 
of an electron between the orbital v on the site 
i, and fi on j, together with the orbital-dependent on- 
site direct- Coulomb (t/y,/i) and the exchange (J,/,^) in- 
teractions. The ratio of the effective interactions (so 
called U representing the orbital average of U^^^) to 
the bandwidth (or averaged transfer t) substantially (ap- 
proximately 50%) increases in the order from LaFePO, 
LaFeAsO, BaFesAsa, FeTe to FeSe^. This is physically 
well understood by the increasing distance h between 
an iron layer and the neighboring pnictogen/chalcogen 
atoms, which alters the chemical bonding of the Fe 3d 
Wannier orbitals with the pnictogen/chalcogen p orbitals 
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from more covalent to more ionic in this order. This 
makes the Fe 3d Wannier orbitals more localized and the 
screening by pnictogen/chalcogen p and other orbitals 
less effective. Therefore, the family dependence offers a 
good example of the interaction control. 

The interaction ?7^,;y strongly depends on the orbitals 
/i, v, and compounds. Nevertheless, all the derived ab 
initio models are reasonably well reproduced from that 
of a particular compound, say LaFeAsO, by a single pa- 
rameter A which scales all the interactions, U^^^ and J^.^ 
uniformly (see Supplemental Materials S.2). In fact, A is 
a measure of the electron correlation. 

Result. — In FigH] (a), we show the ordered magnetic 
moment of the AF order TO(gpoak) (red circles) as a func- 
tion of A, calculated by mVMC with the extrapolation 
to the thermodynamic limit. Here, A = 1 represents the 
ab initio model for LaFeAsO, while by uniformly scaling 
all the interaction strength by A, the interaction control 
is monitored. We also plot the mVMC result of the real 
ab initio models for LaFePO, BaFe2As2, and FeTe, as 
balck- framed yellow diamonds, at corresponding A val- 
ues (see Supplemental Materials S.2)[23|- The result for 
the ab initio models of FeTe shows the double-stripe or- 
der degenerate with the simple stripe with a nearly equal 
ordered moment. All the other results shown here indi- 
cate the stripe-ordered ground states. The agreement 
with the experimental results with the correct AF quan- 
tum critical point proves the accuracy of the present 
scheme. Similar quantitative agreements with the ex- 
periments have been suggested in earlier studies 



(a) 



The uncertainties in these A-scaling plots for results of 
the experiments (blue asterisks) and the ab initio models 
are indicated by horizontal error bars. Here, the uncer- 
tainty arises from the fact that two independent ab initio- 
model derivations give a slight difference (~ 5%) to- 
gether with the error caused by details of orbital varia- 
tions ignored in the uniform A scaling ('^ 5%). 

Quantum critical point of the AF transition appears 
at around (A ^ 0.95). The overall agreement among the 
A-scalcd models, ab initio models and the experiments 
show that the material dependence of the magnetism is 
well described by the variation of the correlation strength 
represented by the single parameter A. 

In addition to the interaction control, they may be sub- 
stituted by other elements so that carriers are doped by 
holes or electrons (namely the filling control) as compared 
to the mother materials at dP. This simultaneous possi- 
bility of the interaction and filling controls makes the 
iron-based superconductors valuable and unique in com- 
parison to the cuprates (mainly, filling control only) or 
organic conductors (mainly, bandwidth control only) [221 ■ 
Furthermore, nearly degenerate five 3d orbitals enable us 
to examine multi-orbital effects, that are absent in the 
cuprates and the organic conductors. 

In the plane of A as well as the doping concentration 
6, the magnetic ordered moment m(gpoak) is plotted in 
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FIG. 1: (color), (a) Magnetic ordered moment m(gpcak) at 
the Bragg wavenumber Qpcak (for the pattern, see right) cal- 
culated by mVMC (red circles). The result at A=l repre- 
sents that of the ab initio model for LaFeAsO while the A 
dependence illustrates the results of the models obtained by 
uniformly scaling the interaction strengths by A. The black- 
framed yellow diamonds indicate the results of the ab initio 
models, where they are plotted at the corresponding A. The 
plotted results are those in the thermodynamic limit after 
the size extrapolation. Experimentally observed values[^-Q] 
are also plotted at corresponding A by asterisks, where FeTe 
shows double-stripe and others show stripe orders in agree- 
ment with our ab initio results, (b) Magnetic ordered mo- 
ment in the plane of A and doping concentration 5. The data 
are plotted in the cross-sections for A=l as well as for 6=0. 
It shows a peak at (5=— 1.0) and decreases monotonically 
over d^ (5=0), which forms a large half-dome structure peaked 
at d'' . The mVMC results of the ab initio models are shown 
by black-framed yellow diamonds. The green and blue shaded 
regions represent the G-type and stripe-type AF orders, re- 
spectively, as is illustrated in the right panel of Fig[T] (a). 
There exist two first-order transitions (black dashed lines), 
one indicated by the jumps in the ordered moment around 
5~0.17 and the other at the transition between the G-type 
and stripe around 5^—0.22, which signals large charge fiuc- 
tuations under phase-separation effects. In the present short- 
ranged-interaction model, the phase separation indeed occurs 
in the gray shaded regions. 



Fig. [T] (b). Here, as in A, i5 monitors the doping ef- 
fect, if the real material could be purely doped without 
changing other parameters such as the transfers. By this 
monitoring, we can elucidate physics of the filling con- 



trol. As we clariiy later, we caution that it does not 
simply mean the substitution effect of the real material, 
because the real substitution may change other parame- 
ters. Although the region near the configuration has 
not been experimentally realized so far, our monitoring 
of the filling control predicts that the iron-based super- 
conductors around 6=0 (d^) are located in the foot of 
a large dome centered at S=—1.0 (d^)- It was formerly 
believed that metallic "valleys" at non-integer fillings in- 
tervene the insulators formed at each integer filling (see 
for instance Fig. 65 of Ref. In marked contrast, the 

monotonic decrease in TO(gpeak) from d^ over d^ supports 
that the electron correlation of the iron-based supercon- 
ductors around d^ emerges just as a proximity effect of 
the d^ Mott insulator without a good metallic region be- 
tween d^ and d^ for the retained interaction strength. 

Indeed the orbital-resolved momentum distribution 
riuik) in Figl2] shows that the clear jump manifested 
at the Fermi surface of good metals becomes less and 
less clear with decreasing electron concentration from d^ 
(5=0) toward d^ {5=—l). The reason for the strong inco- 
herence toward d^ is that all the orbitals become nearly 
half filling at d^ . Such a strong commensurability leads 
to a solid Mott insulator. The large dome structure is 
a consequence of an overwhelming strong proximity of 
the d^ (half- filled) Mott insulator, which blankets the 
d^ commensurability. Furthermore, signatures of several 
first-order transitions seen in FiglT] (b) in the filling con- 
trol signal large charge fluctuations with tendency toward 
the phase separation and incoherence. 

Another important aspect found in the doping depen- 
dence is the switching in the magnetic order; the tran- 
sition from the stripe to G-type orders ai 6 ^ —0.22 
[Fig.l (b)], which is explained by the geometrical frustra- 
tion effects: The ratio of the diagonal next-nearest (i') 
to the nearest-neighbor (i) transfers, i.e., t' /t measures 
the frustration effect. In fact, dyz/dzx and dx^-y^ or- 
bitals have t' /t ~ 1.0, while for dxY and dz^ orbitals 
t'/t - 0.1 0, where the orbitals are labeled in the repre- 
sentation of the folded Brillouin zone Around d^, the 
frustrated dx^-y^ is dominant because of the differenti- 
ated pinning to half filling, while near d^ , all the orbitals 
are pinned to half filling and dxy and dz^ become the 
more important players because of the larger effective 
interactions, which favors the G-type order. Nearby G- 
type order is suggestive of proposed nodes in the pairing 
symmetry of heavily hole doped KFe2As2[2l,[23l. In this 
connection, the double-stripe type order at (7r/2, tt/2) ex- 
perimentally observed for FeTe is also accounted for by 
the involvement of the less frustrated dz^ orbital. 

To further understand the proximity effect of d^ Mott 
physics, we monitor the effects of Coulomb (U) and ex- 
change (J) interactions separately by the scaling parame- 
ters Xu and Xj, respectively. Figure [3] reveals that U and 
J comparably contribute at d^ in enhancing the ordered 
magnetic moment. The Hund's rule coupling J retains 
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FIG. 2: (color). Filling dependence of orbital-resolved mo- 
mentum distribution n^{k) = {c]^ ^c^,k) plotted for the "un- 
folded" Brillouin zone (BZ), where the creation and anni- 
hilation operators of an electron at the orbital v and the 
wavenumber k are denoted by j. and c,^,k, respectively. In 
the bottom left, the unfolded BZ (with the coordinates kx 
and ky) and folded BZ (with kx and ky) are depicted by the 
red and black lines, respectively. In the bottom right, the 
Fermi surfaces of the LDA band structures are shown by the 
red curve, which can be identified in the corresponding sharp 
changes in nv{k) for v = dyz/zx and djjj-2_y2 on the top pan- 
els. Note that the orbitals are represented according to the 
folded BZ. Fade out of sharp boundaries in nv{k), represent- 
ing the Fermi pockets, are seen especially for dyz/zx and 
dx2_Y2 together with d^-i, which signals strong renormal- 
izations of quasiparticles with bad metallic behavior when 5 
decreases progressively to negative. 



the metallic but AF order even far away from half filling 
(d^) contrary to the quick collapse of the AF order in the 
cuprates. The present result conforms with a suggestion 
of an importance of J away from half filling 27| . We see 
below that J between dx'^-y'^ and dz'^ orbitals is crucial 
in the iron families. 

Around d^, large spin and orbital fluctuations are 
prominent with the nearby magnetic quantum critical 
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FIG. 3: (color online), (a) Magnetic ordered moments in the 
plane of Xu and Xj for system size A^s =6x6, which is 
expected to be close to the thermodynamic limit. Here, Xj 
(Xu) scales exchange (Coulomb) interactions from the ab ini- 
tio model for LaFeAsO at Xu = Xj = 1. The sharp change 
in the ordered magnetic moment is triggered by a synergy 
effect of the direct Coulomb (Xu) and the exchange (Aj) in- 
teractions and the two interactions comparably contribute, 
(b) Orbital-resolved double occupation defined by proba- 
bility of simultaneous occupation of up and down spin elec- 
trons on the f orbital at the same site, (c) Orbital occupation 
defined by the averaged density of electrons on the orbital 
u as a function of A at filling S — O(d^). Antiferromagnetic 
phase is represented as shaded region. System sizes 8x8 for 
(b) and (c) are sufficient to well represent these quantities in 
the thermodynamic limit. Near the realistic parameter A = 1, 
strong correlation of the dx2-Y^ orbital synergetically starts 
involving the dz2 orbital through the Hund's rule coupling, 
which drives the into the state close to half filling n — 1 
with reduced D. 



point. Figures |3] (b) and (c) illustrate that the dx^-y^ 
orbital is under the strongest correlation effects indicated 
by the smallest double occupation [panel (b)] and a re- 
markable pinning around half filling [(c)] even for rela- 
tively small A, though the effective Coulomb interaction 
on the dx2-Y2 is the smallest even at A = I 28|, i.e., on- 



site intraorbital interaction for the dx^-y^ orbital is 1.68 
eV in comparison to 2.24-2.75 eV for other intraorbital 
values. (See Table I in Supplemental Materials.) This 
puzzle is solved by the largest bare density of states at 
the Fermi level (sec Fig. 3 of Ref. Q), which sensitively 
allows the electron correlation through the electron-hole 
polarizations. When A increases around A = 1 (real- 
istic for LaFeAsO), the orbital quickly follows up 
the dx2-y2, because the dx'^-y^ orbital has the largest 
Hund's rule coupling J with dz^ (0.43 eV in comparison 
to 0.23 eV with dxy and 0.35 eV with dyz/zx)- The 
moment at the dz^ orbital grows, dragged by dx^-y^ 
through J. We confirmed that such orbital differentia- 
tions occur commonly in the ab initio models for other 
families |15|.The strong crossover from low to high mag- 



netic moments by changing from LaFeAsO, BaFe2As2 
through FeTe is accounted by this interplay of dx^-y^ 
and dz2 orbitals. 

One might suspect that the disappearance of the mag- 
netic order in hole doped Bai_a;Ka;Fe2As2 at a; 0.4[25j 
appears to contradict the present result. However, K 
doping does not work as a simple filling control but yields 
more complicated effects, because the Fe square lattice 
shrinks upon the K doping. This simultaneously en- 
hances the bandwidth (nearly 10%) and works also as the 
bandwidth control driving into weak correlation regime. 
At (5 = — 1 and A = 1, we obtained the d^ Mott insu- 
lator with a large Mott gap 2.2eV) with a high or- 
dered moment (~ 2.5/xb)- It implies a connection with 
a G-type AF insulator in isostructural but d^ compound 
LaMnP0[2^. 

Summary. — Wc have shown that the underlying prox- 
imity effects of the d^ Mott insulator governs the elec- 
tronic structures of iron-based superconductors, which 
are located around d^ filling. When the Mott proximity is 
weakened, strong orbital and spin fluctuations take place 
before the verge to weakly correlated metals. This is in- 
deed the high- Tc superconducting region around the d^ 
filling. It shares a common character with the cuprates, 
while a new aspect here is an involvement of the orbital 
fluctuation and differentiation under the big umbrella of 
the d^ Mott insulator. Although d^ is commensurate 
filling, the Dirac nodes in the band dispersions main- 
tain the metal for unexpectedly wide region (2^, [s^], 
while the interactions develop the antiferromagnetism. 
The AF quantum critical points accompany the strong 
crossover of the orbital polarization as well and the su- 
perconducting mechanism has to be clarified under this 
circumstance. We propose to focus particularly on the 
role of the dx'^-y^ and dz^ orbitals as the leading play- 
ers, in contrast to the dyz/dzx orbitals in the nesting 
picture. In the present results, the d^ dominance washes 
away the subtlety of the nesting. Experimental clarifi- 
cation of the large dome structure in the full bandwidth 
and filling controls in the range from d^ to d^ is highly 
desired to clarify the whole perspective of the electron 
correlation effect in the iron-based superconductors. 
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Supplemental Materials 



s.l DERIVATION OF EFFECTIVE 
LOW-ENERGY MODELS BY AB INITIO 
PROCEDURE 

To understand the normal state of iron-based super- 
conductors, we employ the Multi-energy-scale Ab ini- 
tio scheme for Correlated Electrons (MACE), which is 
a hybrid method of ab initio calculations by combining 
density functional theory with an accurate low-energy 
solver [Sl|. We first compute the global electronic struc- 
ture by the local density approximation and eliminate 
the bands far away from the Fermi level by a downfolding 
procedure after renormalizing the effects of these elimi- 
nated bands. The renormalization consists of 1) the con- 
struction of the Wannier orbitals for the target bands 
and the derivation of the tight-binding Hamiltonian in 
the Wannier basis followed by 2) the computation of 
the effective interaction for the target Wannier orbitals 
screened by the eliminated bands. This partial screen- 
ing is computed by following the standard constrained 
RPA method Sl| . After deriving the effective low-energy 
model, 3) it is solved by solvers such as the mult- variable 
variational Monte Carlo method (ssl . [s3| whose frame- 
work is outlined in S.3. This hybrid method has been 
applied to various materials that exhibit strong electron- 
correlation effects including transition metal compounds 
and organic conductors (for a review, see reffsH). In the 
present case, the target bands consist mainly of the Fe 
3d five orbitals. Detailed derivation of thus obtained 3D 
low-energy effective models is given in refslS 



and [83. 

We further downfold the 3D model into the 2D model 
by computing the scree ning effect from the layers other 
than the target laver(s3.IS3l|. The validity of our effective 
model in 2D is confirmed by the fact that the magnetic 
ordered moment is consistent with the experimental val- 
ues as shown in Fig. 1. 

Our 2D low-energy model is written for ten-fold de- 
generate Fe-3d orbitals in a unit cell containing two Fe 
atoms in the form 



(S.l) 



where c^^^,„ (oi^^^a) is a creation (annihilation) operator 
of an electron with s pin a on the i^th maximally local- 



Ji,i,v,tJ. ^'I'e screened Coulomb and exchange interactions, 
respectively. We use the transfer integrals up to the fifth 
neighbors, which well reproduce the LDA band struc- 
tures. Offsite interactions are ignored since those are less 
than a quarter of the onsite parameters and we take the 
site-independent notation f/j,^^ and Jv,fi instead oiUi^i^^^^ 
and J, 



respectively hereafter. Our results of the ab 
initio dimensional downfolding for each compound all in- 
dicate that the screenings from the other layers reduce all 
the Coulomb interactions uniformly [s3|- Based on our 
calculated ab initio results, we subtract the constant val- 



ues from the U^n in 3D models 



as follows; 0.44 eV 



for LaFcPO, 0.41 eV for LeFcAsO, 0.38 eV for BaFcaAsa, 
and 0.40 eV for FeTe, respectively. We note that the ex- 
change interactions J^^^ are not significantly changed by 
the screenings from the other layers [s3l- All the inter- 
actions parameters that we used in this paper are shown 
in Table S I. 

The double counting of the electron correlation already 
considered in the local density approximation (LDA) is 
carefully removed with adjusting of the orbital levels so 
that the Hartree solution of effective models conforms 
with the LDA solution for real materials [s^. By us- 
ing this method, we estimate the double counting part 
proportional to the Coulomb interactions and confirm 
that the inter-orbital part is much smaller than the intra- 
orbital part. Therefore neglecting the inter-orbital part 
hardly changes the present results. Furthermore, by fol- 
lowing the standard way for estimating the double count- 
ing [S6| , it is confirmed that the double counting of the 
exchange correlation is also negligible |S3{ . To examine 
the doping effects, we change the number of electrons 
without changing the band structures, i.e., we maintain 
the orbital levels at those of the undopcd [d^) case. 



S.2 DISCUSSION ON ADEQUACY OF SINGLE 
PARAMETER SCALING BY A 



Although the effective Coulomb interaction U^.^j. rather 
strongly depends on combinations of electronic orbitals 
V, fj,, the ratio to the 3d bandwidth has a clear system- 
atic and uniform change from compounds to compounds 
with relatively small variance. We introduce the domi- 
nant transfer by f = | I]i.=i,5 ^max,i/, where tmax.iy is the 
largest transfer among pairs between an orbital v at an 
iron site and any 3d orbitals at its nearest-neighbor iron 
site. Then the scaled ratio of the interaction strength 
is defined by the dimensionless parameter U^^ = U^fj,/t. 
A single parameter to characterize the ratio of the in- 
teraction to the bandwidth normalized by the value for 
LaFeAsO is 



ized Wannier orbital |S5i at the i-th. site, ti 



contains 



A = 



single-particle levels and transfer integrals. Uis^v.ii and 



(S.2) 
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Table S I: Effective on-site Coulomb (U) / exchange (J) interactions between two electrons on the same iron site in the ab initio 
model (in eV). 



LaFePO 






U 












,7 








XY 


YZ 


Z^ 


ZX 


X^ - Y^ 




XY 


YZ 




ZX 


X'-^ - Y^ 


XY 


2.54 


1.36 


1.34 


1.36 


1.33 


XY 




0.45 


0.54 


0.45 


0.20 


YZ 


1.36 


1.98 


1.53 


1.20 


1.02 


YZ 


0.45 




0.32 


0.36 


0.31 




1.34 


1.53 


2.37 


1.53 


1.02 


Z^ 


0.54 


0.32 




0.32 


0.37 


zx 


1.36 


1.20 


1.53 


1.98 


1.02 


ZX 


0.45 


0.36 


0.32 




0.31 


X^ - Y^ 


1.33 


1.02 


1.02 


1.02 


1.24 


x^ - y'^ 


0.20 


0.31 


0.37 


0.31 




LaFeAsO 






U 












,7 








XY 


YZ 


Z^ 


ZX 






XY 


YZ 




ZX 


X^ - Y^ 


XY 


2.62 


1.39 


1.37 


1.39 


1.50 


XY 




0.46 


0.57 


0.46 


0.23 


YZ 


1.39 


2.02 


1.56 


1.21 


1.11 


YZ 


0.46 




0.33 


0.37 


0.35 




1.37 


1.56 


2.43 


1.56 


1.10 


Z^ 


0.57 


0.33 




0.33 


0.42 


ZX 


1.39 


1.21 


1.56 


2.02 


1.11 


ZX 


0.46 


0.37 


0.33 




0.35 


X^ - Y^ 


1.50 


1.11 


1.10 


1.11 


1.50 


x^ - y'^ 


0.23 


0.35 


0.42 


0.35 




BaFo2As2 






U 












,/ 








XY 


YZ 


Z^ 


ZX 


X^ - Y^ 




XY 


YZ 


Z^ 


ZX 


X^ - Y^ 


XY 


2.80 


1.56 


1.61 


1.56 


1.78 


XY 




0.48 


0.60 


0.48 


0.26 


YZ 


1.56 


2.26 


1.83 


1.39 


1.34 


YZ 


0.48 




0.36 


0.40 


0.41 


Z^ 


1.61 


2.83 


2.90 


1.83 


1.39 


Z" 


0.60 


0.36 




0.36 


0.50 


ZX 


1.56 


1.39 


1.83 


2.26 


1.34 


ZX 


0.48 


0.40 


0.36 




0.41 


X^ - Y^ 


1.78 


1.34 


1.39 


1.34 


1.91 


x^ - y'^ 


0.26 


0.41 


0.50 


0.41 





FeTe 






U 












J 








XY 


YZ 


Z^ 


ZX 


X^ - Y^ 




XY 


YZ 


Z^ 


ZX 


X'-' - Y^ 


XY 


3.44 


1.94 


2.10 


1.94 


2.64 


XY 




0.49 


0.68 


0.49 


0.34 


YZ 


1.94 


2.48 


2.16 


1.63 


1.89 


YZ 


0.49 




0.37 


0.37 


0.49 


Z^ 


2.10 


2.16 


3.44 


2.17 


2.04 


Z" 


0.68 


0.37 




0.37 


0.66 


ZX 


1.94 


1.63 


2.17 


2.48 


1.89 


ZX 


0.49 


0.37 


0.37 




0.49 


X-' -Y^ 


2.64 


1.89 


2.04 


1.89 


3.19 


X^ - Y'^ 


0.34 


0.49 


0.66 


0.49 





where C/lfaoi/^' ^^"^ denominator is the orbital depen- 
dent C/i/^ for LaFeAsO. Then A monitors the orbital av- 
eraged "U/t" scaled by the value for LaFeAsO. We note 
that the values of A for each materials are slightly differ- 
ent from those of ref. |S3| . because In ref. [ssl . we have per- 
formed the dimensional downfolding only for LaFeAsO 
and have assumed the reduction of Coulomb interactions 
is independent of materials. In this paper, we have per- 
formed the dimensional downfolding for each material, 
not only for LaFeAsO, but also for LaFePO, BaFe2As2 
and FeTe to derive the ab initio two-dimensional mod- 
els, and we have precisely determined the reductions of 
Coulomb interactions depending on materials. The re- 
ductions slightly depend on materials as one sees in Sec- 
tion S.l. Small material dependences of the reduction 
(within 10% of the whole reduction) is the reason why the 
estimation of A is slightly different from that of ref. |S3. 
We note that this slight difference did not essentially al- 
ter out results. 

Since the band structure of the Fe 3d bands near the 
Fermi level is similar for all the iron-based superconduc- 
tors with a similar bandwidth around 4.5 eV, the ratio A 
is the only prominent parameter to identify the material 
character of each compound and it is indeed a measure 
for the electron correlation strength. 

To see whether the single parameter A well repre- 



sents the material dependence, we introduce quanti- 
ties luiyfi = C/iyp/C/LFAOi^p, lu = ^J2u,ti^'^f^/^^FAOiJfi 
Iji^fj, = Jvii/ JhFAOvin and Ij = Jvfj./ JhFAO^tJ. 

with Ji/^ = Jv^/i. If /[//A and /j/A are close to unity, it 
supports that the interactions of each material is well 
described by scaling the interactions of the reference 
LaFeAsO. Further support can be obtained by examining 
the variances 



(S.3) 
(S.4) 



In Fig. ^ Ijj, and Ij for several different materials at 
the corresponding A indeed show that /[//A and Ij/X 
are close to unity for all the materials. This means that 
interaction parameters are well scaled uniformly with the 
single parameter A for different families. A relatively 
large deviation from I = X with a substantial variance 
is seen for Ij of FeTe. However, the magnetic moment 
of FeTe is more or less already saturated because of the 
large J7, and this variance does not seriously influence 
the prediction on the moment. 



=D 1 



FeTe 



BaFe2As2 




+ :LaFeAsO 



1.5 



FeTe 



BaFe2As2 




Fig. S 1: Scaling plot of interaction parameters for several 
different families. Definition of / and A are given in the text 
and the error bars show A/y and A/j, the variance (standard 
deviation) of lu and Ij, respectively. The single parameter 
scaling by A works well if the data are on the line I = X with 
small variances. We admit roughly 5% error in the assign- 
ments of the real materials. 



S.3 MULTI- VARIABLE VARIATIONAL MONTE 
CARLO METHOD 

Our variational wave function is defined as 

IV) =7'G7'j/:^=°|0pair), (S.5) 

where Vq and Vj are the Gutzwiller and Jastrow fac- 
tors, respectively [s^j- The Gutzwiller factor punishes 
the double occupation of electrons on the same orbital 
defined by the orbital-dependent variational parameters 
5f as 



Vg = exp - ^ g^rii^^rii^i 



where rii^a- is the density operator for the orbital v with 
the spin a at the i-th site, and the Jastrow factor takes 






J7 























Fig. S 2: The translational symmetry employed for tiling by 
1x4 supercell structure. The labels 1-4 denote the sublattice 
indices for the 1x4 supercell structure for /ij . The arrows in- 
dicate the pattern of the staggered magnetic ordered moment 
realized in the present result of the double-stripe structure. 



into account the intcrorbital and/or intersitc charge cor- 
relations by the site- and orbital-dependent variational 
parameters vijyu as 



Vi = exp 



(S.6) 



In the present study, we take into account up to the near- 
est neighbor for the intersitc pair including the on- 
site intcrorbital pairs by considering all the combinations 
of orbital dependent correlations. 

The spin quantum- number projection C^^^ restores 
the SU(2) spin-rotational symmetry and generates a state 
with the correct total spin 5=0. The one-body part 
0pair) is the generalized pairing wave function defined as 



N/2 



|0) with fij being the vari- 
ational parameters. Here, Ns and N are the system size 
(number of sites) and the number of electrons, respec- 
tively. In this study, we allow fij to have 2x2 supercell 
structure or equivalently we have 2 x 2 x 5^ x iVg indepen- 
dent variational parameters. All the variational parame- 
ters are simultaneously optimized by using the stochastic 
reconfiguration method |S7l.lS8{. The variational function 
I'ip) in Ea. (|S.5|) can flexibly describe paramagnetic metals 
and the stripe type AF and G-type AF phases, as well 
as superconducting phases and correlated metals. In the 
present calculation we also examined the possibility of 
the double stripe type order realized for FeTe by modi- 
fying fij with the 1x4 supercell (sec Fig. S[5]), while we 
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have restricted to the normal state for the metals. For 
the accuracy of the present method see refls?!. 

The magnetic ordered moment is defined by m{q)^ = 
3% E^,J{S^SJ)e''^(^^--^^ for Ar,-site iNj2 unit cell) sys- 
tem with the periodic boundary condition and q is the 
momentum at the peak, which is gpoak = (OiTI") ^md 
(tt, tt) and (7r/2,7r/2) for the stripe, G-type and double- 
stripe orders, respectively. In this definition, the sat- 
urated magnetic moment for the classical Neel state is 
given hy m — 4 fiB- The double occupation is defined 
by = -^^iini.iy^^rii^jj^i) and the orbital dependent 
momentum distribution n^(k) and orbital-dependent fill- 
ing are defined by = ^^X^fc^-fC^) ^^'^ n^{k) = 
X]cr('^t cr(^)'^!^.ff(^))- the physical quantities are pre- 
sented after proper size extrapolations to the thermody- 
namic limit unless explicitly stated otherwise. 
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